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We present an overview of a modeling environment, consisting of a coupled atmosphere-wildfire 
model, utilities for visualization, data processing, and diagnostics, open source software repositories, 
and a community wiki. The fire model, called SFIRE, is based on a fire-spread model, implemented 
by the level-set method, and it is coupled with the Weather Research Forecasting (WRF) model. A 
version with a subset of the features is distributed with WRF 3.3 as WRF-Fire. In each time step, 
the fire module takes the wind as input and returns the latent and sensible heat fluxes. The software 
architecture uses WRF parallel infrastructure for massively parallel computing. Recent features of the 
code include interpolation from an ideal logarithmic wind profile for nonhomogeneous fuels and ignition 
from a fire perimeter with an atmosphere and fire spin-up. Real runs use online sources for fuel maps, 
fine-scale topography, and meteorological data, and can run faster than real time. Visualization pathways 
allow generating images and animations in many packages, including VisTrails, VAPOR, MayaVi, and 
Paraview, as well as output to Google Earth. The environment is available from openwfm.org. New 
diagnostic variables were added to the code recently, including a new kind of fireline intensity, which 
takes into account also the speed of burning, unlike Byram's fireline intensity. 

1. INTRODUCTION 

A model that can be used by others than the authors requires an extensive infrastructure far beyond the simulation 
code itself. Social aspects are just as important as technical ones, if the model is to be used by humans. For this 
reason, the OpenWFM.org environment consists of multiple components, many under continued development: 

• The wildland fire simulation code SFIRE coupled with a mesoscale atmospheric simulation code, the Weather 
Research and Forecasting (WRF) model (Sec. 2). 

• Facilities for acquisition of data from online databases, and preprocessing of fine-resolution data (described 
separately in Beezley et al. (2011); see also Beezley (2011)). 

• Documentation and diagnostic utilities (Sec. 3) 

• Visualization, clients, web portals, and social networks (Sec. 4) 

* Paper 6.4, Ninth Symposium on Fire and Forest Meteorology, Palm Springs, CA, American Meteorological Society, October 201 1 . Available 

at http : //ams . confex . com/ams/9FiRE/webprogr am/Paper 1922 7 7 .html. This research was supported by NSF grant AGS-0835579 
and by NIST Fire Research Grants Program grant 60NANB7D6144. 
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Figure 1 : An overview of the wildfire simulation environment (Mandel et al. 201 0). 



2. COUPLED FIRE-ATMOSPHERE MODEL 
2.1 . WRF coupled with fire spread model 

Wildland fire is a complicated multiscale process. Fortunately, a practically important range of wildland fire behavior 
can be captured by the coupling of a mesoscale weather model with a simple 2D fire spread model (Clark et al. 
1 996a, b, 2004; Coen 2005). Weather has a major influence on wildfire behavior; in particular, wind plays a dominant 
role in the fire spread. Conversely, the fire influences the weather through the heat and vapor fluxes from burning 
hydrocarbons and evaporation. The buoyancy created by the heat from the fire can cause tornadic strength winds, 
and the wind and the moisture from the fire affect the atmosphere also away from the fire. It is well known that a 
large fire "creates its own weather." 

The code combines the Weather Research and Forecasting Model (WRF) with the semi-empirical fire spread 
model (Rothermel 1972), with the fire spread implemented by the level-set method. Alternatives to the Rothermel's 
model are in progress. The model is able to run faster than real time on several hundred cores, with the fire model 
resolution of few meters and horizontal atmospheric resolution on the order of 100m, for a large real fire (Jordanov 
et al. 2011). 

Kochanski et al. (2010) compared simulation results with measurements on the FireFlux grass fire experiment 
(Clements et al. 2007). Dobrinkova et al. (201 1 ) simulated a fire in Bulgarian mountains using real meteorological 
and geographical data, and ideal fuel data. Beezley et al. (2010) simulated afire in Colorado mountains using real 
data from online sources. Jordanov et al. (2011) simulated a large fire in Bulgaria with real data, including a fuel 
map derived from satellite mesurements, and compared the results with the actual fire area. 

Further details, references, and acknowledgements can be found in Mandel et al. (2009, 201 1b, a). See also 
Fig. 1 for an overview poster. 
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2.2. Related work 

Wildland fire models range from tools based based on fire spread rate formulas (Rothermel 1972, 1983), such as 
BehavePlus (Andrews 2007) and FARSITE (Finney 1998), suitable for operational forecasting, to sophisticated 3D 
computational fluid dynamics and combustion simulations suitable for research and reanalysis, such as FIRETEC 
(Linn et al. 2002) and WFDS (Mell et al. 2007). BehavePlus, the PC-based successor of the calculator-based 
BEHAVE, determines the fire spread rate at a single point from fuel and environmental data; FARSITE uses the fire 
spread rate to provide a 2D simulation on a PC; while FIRETEC and WFDS are physically more accurate and run 
much slower than real time. 

The level set-method was used for a surface fire spread model in Mallet et al. (2009). Filippi et al. (2009) coupled 
the atmospheric model, Meso-nh, with fire propagation by tracers. Tiger (Mazzoleni and Giannino 2010) uses a 2D 
combustion model based on reaction-convection-diffusion equations and a convection model to emulate the effect 
of the fire on the wind. FIRESTAR (Morvan and Dupuy 2004) is a physically accurate wildland fire model in two 
dimensions, one horizontal and one vertical. UU LES-Fire (Sun et al. 2009) couples the University of Utah's Large 
Eddy Simulator with the tracer-based code from CAWFE. See the survey by Sullivan (2009) for a number of other 
models. 



2.3. New features 

Recent features of the code include interpolation from an ideal logarithmic wind profile for nonhomogeneous fuels 
(Mandel et al. 2011a), important for simulation at a variety of scales, and ignition from a fire perimeter with an 
atmosphere and fire spin-up (Kondratenko et al. 201 1). The perimeter ignition is important for practical application 
and the technique will be adapted also for data asimilation, where the atmosphere state needs to be adjusted when 
the fire state changes due to new data. 

New features which were added to the code since the paper Mandel et al. (201 1 a) include computation of f ireline 
intensities. Byram's f ireline intensity (Byram 1959) is the heat produced per unit length of the fireline in unit time 
(J/m/s) in the so-called flaming zone behind the fireline. Hence, it is given by 

I = HRw (1) 

where H (J/kg) is the heat contents of the fuel, R (m/s) is the spread rate, and w (kg/m 2 ) is the fuel amount that 
burns in the flaming zone. In practice, the fuel amount burned to is estimated as a fixed fraction of the fuel load to 
(kg/m 2 ). 

Though Byram's fireline intensity is routinely used for practical guidance, it does not depend on the speed of 
burning. However, if the fuel burns slowly, much of the burning takes place at a distance from the fireline and thus 
may not contribute much to the severity of the fire, while a fast burning fuel will release its heat close to the fireline. 

For this reason, we introduce a new concept of fireline intensity as the amount of heat generated by the 
advancing fireline from the newly burning fuel only, in a small unit of time. Assume that the fuel fraction after 
ignition decreases exponentially with the time t from ignition, as 

where T { is the fuel burn time, i.e., the time when 1 - e _1 w 64% of fuel has burned. Straightforward application of 
calculus then shows that the new fireline intensity is 

J=^ (J/m/s 2 ). (2) 

Unlike Byram's fireline intensity (1), the new fireline intensity, given by (2), takes into account the effect that a 
faster burning fuel will create a more intense heat concentrated at the fireline. The reason why the time unit is 
squared is that over a longer time, the fireline advances over a longer distance, and the fuel it advances over also 
has longer to burn. 

The code computes the fireline intensities and the reaction intensity (which is the same as the released heat flux 
intensity, J/s/m 2 ) for the simulated fire to estimate its severity. The fireline intensities are computed from the fire rate 
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Welcome to the WRF-Fire pages at openwfm.org, the open wildland fire modeling community. WRF-Fire is a freely available wildfire modeling code, which combines a standard weather 
modeling code with a surface fire model. 

In the pages listed below, you can find support information, documentation, and quick guides, called Howtos, for WRF-Fire, as well as links to the download of the latest version. If you 
have any questions or if there is something you would like implemented in WRF-Fire, please let us know by email to jan.mandel@gmail.com 
Main page: WRF-Fire 

Pages in category "WRF-Fire" 

The following 35 pages are in this category, out of 35 total. 

c 

■ Changes in WRF-Fire 3.3 release 



■ WRF-Fire development notes 

■ WRF-Fire documentation 

H 

■ How to convert data for Geogrid 

■ How to diagnose fuel properties in WRF-Fire 

■ How to get WRF-Fire 

■ How to interpret WRF variables 

■ How to load WRF files in Matlab 

■ How to run WRF-Fire 

■ How to run WRF-Fire with real data 

■ How to run the standalone fire model in WRF-Fire 



H cont. 

■ How to visualize WRF-Fire output in Google Earth 

■ How to visualize WRF-Fire output in Matlab 

■ How to visualize WRF-Fire output in Mayavi2 

■ How to visualize WRF-Fire output in VAPOR 

■ How to visualize WRF-Fire output in VisT rails 

■ How to visualize vertical profiles from WRF in Matlab 



■ List of SFIRE pages 

■ List of WRF-Fire pages 



S cont. 

■ SFIRE variables 

■ WRF-Fire support 



■ Namelist.fire 

■ Name list, input 



i WRF-Fire publications 



■ Talk:WRF-Fire development notes 

■ Template:WRF-Fire-commi1 

■ Template:WRF-Fire-file 

U 

■ Coupled WRF and SFIRE model user's guide 

V 

■ Vertical wind interpolation 

■ Visualization in Google Earth 

w 

■ WPS with GeoTIFF support 
WRF-Fire 

■ WRF-Fire ignition 

■ WRF-Fire wish list 



Figure 2: Support and documentation pages on the OpenWFM wiki at http : //www. openwfm. org/wiki/List 
of_SFIRE_pages 



of spread R. Since R is well-defined on the fireline only, the fireline intensities are defined on mesh nodes next to 
the fireline only as well. 

Separate computations are made as a component of a fire danger rating to estimate the severity of a potential 
fire, that is, to answer the question: if a fire breaks out, how bad would it be? These quantities are computed from 
the maximal rate of spread at any location for the wind speed and the slope at that location, and they can be used 
to plot potential fire severity maps. 

3. OPENWFM WIKI 

Fig. 2 shows the list of all wiki pages for the support of the simulation code. The pages starting with "How to" are 
short guides that allow the user do basic things simply. Users can then investigate variations at their leisure. 

Diagnostic utilities, described in the wiki and available from source code repositories, allow the user, e.g., 
generate graphs of fire spread rate as a function of wind speed to aid in debugging fire properties, or draw wind 
profiles at a specified point (the wind speed as a function of height). 

Since most the user's guide (appendix in Wang et al. (2010)) was extracted from the wiki articles anyway, and 
the wiki has much more information, such as the diagnostic utilities, the existing wiki articles were now organized 
into a new, continously updated User's Guide at http : / /www. openwfm. or g/wiki /user s_guide (Fig. 3). This 
required adding only a minimal amount of material. 
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Coupled WRF and SFIRE model user's guide 

(FteSreclea Iron Usees guide) 

This guide applies to the code as available from OpenWFM.org and partly also to WRF-Firc. which is a subset available in WRF release This guide evolves with the 
software and it complements the technical description in J. Mandol, J. D. Bcozlcy. and A. K. Kochanski, Coupled atmosphcrc-wUdland fire modeling with WRF 3.3 and 
SFIRE 2Q1 1. Geo^GientfTtG Model Development (GMD) 4. 591-610. 2011. da-i.10.5134fgmd-4*591-2011®, 
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Figure 4: Visualization in Google Maps web portal, 2009 Harmanli fire, Bulgaria. Simulation data from Jordanov 
et al. (2011). 

The OpenWFM wiki provides a social environment, where users can edit and comment on each other's work. It 
is using the WikiMedia software. 

4. VISUALIZATION, CLIENTS, WEB PORTALS, AND SOCIAL NETWORKS 
4.1 . Google Maps and Google Earth 

We have proposed Google Earth visualization already in Douglas et al. (2006). Since then, Google API has become 
a de-facto standard in wildland fire visualization. Both Google Earth and Google Maps use the same file format, 
KML (Keyhole Markup Language). We use the compressed format, KMZ. A KMZ file is a zipped directory with 
XML code containing drawing instructions and metadata, such as geolocation, and (optionally) image files to be 
displayed on the Earth surface. We use the image files to display the heat flux (Figs. 4 and 5) or fire area (Fig. 6). 
We render the file in Google Maps by posting it on the web and providing its URL to the Google Maps Javascript 
API, in a web page (Fig. 4), or loading the file into the Google Earth application. Google Earth can animate a 
KML file with multiple images as frames. Animation in Google Maps requires custom Javascript coding and is under 
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Figure 5: Visualization in Google Earth client, 2010 Mt Carmel Fire, Israel. 

development. Also, surface images to be visualized in Google Earth and Maps API use a different map projection 
than WRF, so the mesh nodes in WRF and the pixels in the image do not line up, and the simulation values need to 
be interpolated. 

Future extensions include a web portal or a client application based on Google Maps, which can initiate a fire 
simulation on a cluster. 

4.2. VisTrails and Crowd Labs 

VisTrails (Freire et al. (201 1), http : / /www. vist rails . org) is a scientific development environment for workflow 
and provenance management. It transparently stores user interactions within the system to ensure reproducibility of 
the produced data sets and visualizations. VisTrails provides a graphical programming interface that was specifically 
tailored for the generation of custom scientific visualizations. Through a drag-and-drop interface, users may 
interconnect modules that are responsible for reading a file, parsing the data, and performing different visualization 
tasks. VisMashups (Santos et al. 2009) is a framework for the automatic development of simplified user interfaces 
from annotated workflows developed within the VisTrails environment. 

Together, VisTrails and VisMashups, provide the tools to generate customizable visualizations for coupled WRF 
and fire simulations. Demonstrated in Fig. 7 is a fully interactive 3D visualization environment of the wildland fire 
simulation results, running locally on a user's machine. These visualizations support exploration of the wildland fire 
through an interactive camera, toggles for the various ground and atmospheric variables, controls for isosurface and 
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Figure 6: Visualization in Google Earth client, 2007 Witch fire, CA. 



timing components, as well as interactive vector streamlines. 

Providing web services, we improve the accessibility of the coupled WRF and fire simulation results for 
educational purposes and to improve communication and data sharing within the research community. We have 
setup the wildland fire simulation user group within CrowdLabs (Mates et al. (2011), http : //www. crowdiabs . 
org). CrowdLabs is a social web site that provides a scalable environment for collaborative analysis involving 
simulation scientists, visualization experts, and other users alike. Through a web portal, users share and discuss 
VisTrails workflows and VisMashup interfaces to improve their efficacy. Fig. 8 illustrates the online interface for 
CrowdLabs, as well as a shared fire visualization on the site. It remains future work to provide 3D interactions within 
the web-based visualizations, as well as develop the controls to manage simulation executions from the site. 
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